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Abstract. X-ray absorption and emission features arising from the inner-shell transitions in iron 
are of practical importance in astrophysics due to the Fe cosmic abundance and to the absence of 
traits from other elements in the nearby spectrum. As a result, the strengths and energies of such 
features can constrain the ionization stage, elemental abundance, and column density of the gas 
in the vicinity of the exotic cosmic objects, e.g. active galactic nuclei (AGN) and galactic black 
hole candidates. Although the observational technology in X-ray astronomy is still evolving and 
currently lacks high spectroscopic resolution, the astrophysical models have been based on atomic 
calculations that predict a sudden and high step-like increase of the cross section at the K-shell 
threshold (see for instance Ref. [1]). 

New Breit-Pauli R-matrix calculations of the photoionization cross section of the ground states 
of Fe XVII in the region near the K threshold are presented. They strongly support the view that the 
previously assumed sharp edge behaviour is not correct. The latter has been caused by the neglect 
of spectator Auger channels in the decay of the resonances converging to the K threshold. These 
decay channels include the dominant KLL channels and give rise to constant widths (independent 
of n). As a consequence, these series display damped Lorentzian components that rapidly blend to 
impose continuity at threshold, thus reformatting the previously held picture of the edge. Apparent 
broadened iron edges detected in the spectra of AGN and galactic black hole candidates seem 
[2, 3] to indicate that these quantum effects may be at least partially responsible for the observed 
broadening. 


INTRODUCTION 

X-ray absorption and emission features due to inner shell transitions in iron are of prac- 
tical importance in astronomy, owing to the cosmic abundance of iron and to the absence 
of features from other elements nearby in the spectrum. As a result, the strengths and 
energies of iron K shell features can constrain the ionization state, elemental abundance, 
and column density of gas in the vicinity of exotic objects such as active galactic nuclei 
(AGN) and galactic black hole candidates. Although the observational technology in X- 
ray astronomy is still evolving, many such features have been observed. Emission lines 
have been used to constrain the geometrical distribution of gas in compact binary stars 
[4], and to detect black holes in AGN [5]. Iron absorption features have been detected 
from AGN [6], from compact binaries [7], and from galactic black hole candidates [8, 3], 
These divide naturally into bound-bound and bound-free transitions, and much of the in- 
terpretation of the bound-free absorption features has been based on atomic calculations 
such as those of Vemer and Yakovlev [9] or Berrington and Ballance [1] which predict 



a sudden increase in the cross section at the binding energy of the K-shell electron. 

When a photon is sufficiently energetic to promote a K-shell electron to a higher 
Rydberg state, the resulting excited state can decay via Auger channels. In the case of 
Fe XVII, we have these possible decay pathways: 

hv + 2p 6 -» ls _1 2p 6 «p 

f 2p 5 + e~ 1 
\ 2s _1 2p 6 + e~ J 

{ 2p 4 np + e~ > 

2s -1 2p 5 np + e~ > 

2s~ 2 2p 6 np + e~ J 

Eq. (1) represents the participator (KLn) channels where the np Rydberg electron is 
involved in the decay while it is not in the spectator (KLL) channels (Eq. (2)). R-matrix 
calculations have been carried out to analyze the role of the spectator channels. 

In the following sections, we will detail the methodology and discuss the results. 
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THEORY 

The Breit-Pauli R-matrix method (BPRM) is widely used in electron-ion scattering and 
in radiative bound— bound and bound— free calculations. It is based of the close-coupling 
approximation [10] whereby the wavefunctions for states of an A/ -electron target and a 
colliding electron with total agular momentum and parity Jn are expanded in terms of 
the target eigenfunctions 


(3) 

i r j 

The functions Xi are vector coupled products of the target eigenfunctions and the angular 
part of the incident-electron functions, F)(r) are the radial part of the latter and 
is an antisymmetrization operator. The functions are bound-type functions of the 
total system constructed with target orbitals; they are introduced to compensate for 
orthogonality conditions imposed on the f](r) and to improve short-range correlations. 
The Kohn variational gives rise to a set of coupled integro-differential equations that 
are solved by R-matrix techniques [11, 12] within a box of radius, say, r < a. In the 
asymptotic region (r > a), solutions are found using the MQDT method [15, 16]. Breit- 
Pauli relativistic corrections have been introduced in the R-matrix suite by Scott and 
Burke [13], Scott and Taylor [14]. Inter-channel coupling is equivalent to configuration 
interaction in the atomic structure context, and presents a formal and unified approach 
to study the decay properties of both bound states and resonances. 



STRATEGY 


In order to display the different effects of these two types of Auger decays, two calcula- 
tions of the photoionization cross section have been earned out (hereafter referred to as 
Calculations A and B ). In Calculation A, the BPRM method have been used to solve the 
e~ + Fe 17+ system. The following target configurations were considered: 2p 5 , 2s _1 2p 6 , 
and ls~*2p 6 . Only the participator Auger channels described in Eq.(l) were therefore 
included. The orbitals of the target ion were optimized using the code SUPERSTRUC- 
TURE [17]. As it was obviously impossible to include all the 2p 4 np, 2s 1 2p 5 np, and 
2s~ 2 2p 6 np target configurations in our R-matrix model. Calculation B was earned out 
in two steps. Firstly, an R-matrix calculation was carried out adding the contributions 
of KLL Auger decays of the ls~'2p 6 3p resonances (Eq.(2)>. The next step consisted 
in fitting Lorentzian profiles to these 3p resonances in order to obtain the widths and 
areas. These symmetric profiles are due to the dominance of the KLL channels over 
the KLn ones. Actually, these spectator Auger channels produce continuum states that 
cannot be reached by direct photoionization of the ground state causing the Fano q pa- 
rameter to go to infinity (see Nayandin et al. [18] and references therein). As the Auger 
rates of the spectator channels are constant along the Rydberg series, the cross section 
was extrapolated for the higher members of the series using Lorentzian profiles with a 
constant width and areas decreasing as the third power of the effective quantum number, 
n*. The resonance positions were deduced by the usual Ritz formula using the data of 
Calculation A. 


RESULTS AND DISCUSSION 
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FIGURE 1 . Total photoionization cross section in Mb versus photon energy in keV. Solid and dash lines 
are Calculations B and A iespectively (see the text).The sharp K-edge near 7.7 keV in Calculation A has 
disappeared in Calculation B. 




In Fig. 1, a comparison between Calculation A and B is shown. The sharp K-edge 
near 7.7 keV has disappeared, being filled up by the highest members of the Rydberg 
series which are blended together. Note also the differences in the resonance shapes 
and intensities. Because of the constant width, the intensities decrease as l/(n ) in 
Calculation B. In Table 1, a comparison between widths obtained by two R-matrix 
calculations (with and without the spectator Auger channels) is presented for the 3p and 
4p resonances. One can see the significant effect of the spectator (KLL) Auger decays 
on the widths. As the KLn widths (T KLn ) are going down as (n*)~ 3 and the KLL ones 
are independent of n* , the resulting widths (T +kll^ 5 e come constant for the highest 
members of the Rydberg series. Actually, see Table 1, we can conclude that they are 
almost constant for the first few members. 

TABLE 1. Comparison between widths obtained without 
(r KLn ) and with the KLL (spectator) Auger decays (r KLn+KLL ) 
using the BPRM method for the 3p and 4p resonances. 


Resonance 

E (keV) 

F KLn (eV) 

pKL/i+KLL (eY) 

ls -1 2p 6 3p 3 P| 

7.1865 

3.802x 10~ 2 

6.495 x 10 -1 
6.335x 10 -1 

ls-^p^p'P? 

7.1923 

2.165x 10~ 2 

ls-'2p 6 4p 3 P? 

7.4224 

1.386xl(T 2 

6.146x 10 _1 

ls- , 2p 6 4p 1 P‘j’ 

7.4246 

7.952x 10 -3 

6.082 x IQ- 1 


Fig. 2 depicts the cross section obtained in Calculation B convoluted with an in- 
strumental (Gaussian) profile of 10 eV. The ‘edge’ is now shifted down to 7.6 keV. 
The 9p resonances are the highest resolved member of the Rydberg series. Unfortu- 
nately, we have found no experimental data in Fe xvn to compare with our calcu- 
lations. Nevertheless, Farhat et al. [19] have measured the total electron yield of the 
hv + 3p 6 — t 2p -1 3p 6 + e~ photoionization process in neutral argon. The 2p _1 3p 6 n^ 
resonances can decay via dominant spectator Auger channels to the 3p*nl, 3p _1 3p 5 nC, 
and 3s~ 2 3p 6 n£ Ar II configurations. Fig. 3 is the reprint of Fig. 2 in Farhat et al. [19] 
which presents the total electron yield versus photon energy (the resolution is 160 meV). 
The 2p 3/2 threshold is indicated at 248.6 eV but no sharp edge is seen at this energy. Note 
the similar behaviour near threshold with Fig. 2. Actually, Gorczyca and Robicheaux 
[20] used an optical potential within the MQDT method to mimic the effect of the spec- 
tator Auger channels on the L-shell photoionization cross section of neutral argon. This 
technique has been also compared to measurements of the total electron yield near the 
K-shell threshold in neutral neon [21] and oxygen [22]. 

The calculations presented here show that the K-shell photoabsorption cross section 
for Fe XVII does not exhibit a sharp edge at the the binding energy of the K-shell 
electron, but it is shifted and broadened towards lower energies. It is worth pointing 
out that apparently broadened iron edges have been detected in the spectra of AGN and 
galactic black hole candidates [2, 3], and the same physical processes we identify here 
may be at least partially responsible for the observed broadening. 




FIGURE 2. Total cross section obtained in Calculation B (see the text) convoluted with a Gaussian 
profile with a width of 10 eV. 



FIGURE 3. Total electron yield versus photon energy in eV as measured in Ar I by Farhat et al. [ 19]. 
The energy scan was taken over the 2p~W resonances. 


CONCLUSION 

As an important conclusion, this work shows that the sharp edge structures at thresholds 
disappear when the Rydberg series before threshold are affected by spectator Auger 
decays. These inner-shell Auger channels cause constant-width Lorentzian shape reso- 
nances that blend near the threshold changing drastically the structure of the edge. Al- 
though this effect has been measured for the L-edge in neutral argon and for the K-shell 
in neutral neon and oxygen, experimental data are needed for the K-edge in Fe xvii. 
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